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We study the effect of varying strength, 5, of bond randomness on the phase transition of the 
three-dimensional Potts model for large q. The cooperative behavior of the system is determined by 
large correlated domains in which the spins points into the same direction. These domains have a 
finite extent in the disordered phase. In the ordered phase there is a percolating cluster of correlated 
spins. For a sufflciently large disorder 5 > St this percolating cluster coexists with a percolating 
cluster of non-correlated spins. Such a co-existence is only possible in more than two dimensions. 
We argue and check numerically that St is the tricritical disorder, which separates the first- and 
second-order transition regimes. The tricritical exponents are estimated as f3t/i^t ~ 0.10(2) and 
ft = 0.67(4). We claim these exponents are q independent, for sufficiently large q. In the second- 
order transition regime the critical exponents fit/i't ~ 0.60(2) and ut — 0.73(1) are independent of 
the strength of disorder. 



I. INTRODUCTION 

The effect of bond randomness on the critical behav- 
ior of ferromagnetic models is well understood if the 
non-random system has a second-order phase transitior*i. 
Much less is known, however, if this transition is 
discontinuous^. In two dimensions (2d) rigorous results 
asserts that for any type of continuous disorder the tran- 
sition softens into a second order ono^. Recent numeri- 
cal studies of the 2d (/-state random bond Potts model^ 
(RBPM) have shown that the magnetization exponent of 
this model is q-dependent^ and saturatesi*^ for largc-g 
at a possibly exactly known valua^tiS. On the other hand 
the energy exponent, is found to show only a very weak 
variation with q. 

In real 3d systems the effect of bond randomness is 
more complex and here we are lacking rigorous results. 
It is demonstrated experimentally that the isotropic to 
ncmatic transition of nCB liquid crystal turns to second 
order for sufficiently strong disordci*ii. The same type of 
softening effect is found in Monte Carlo simulations for 
the g = 3 and g = 4 Potts models, both for sit o^^d^ and 
bond dilutioni4. With very large computational effort it 
was possible to locate the second order transition point 
and to estimate the critical exponents, which are found q 
dependent both for the magnetization and for the energy 
density. 

The first-order transition regime (for weak disorder) 
and the second-order transition regime (for strong disor- 
der) are separated by a tricritical point the properties of 
which are conjectured to be related to the critical point 
of the random field Ising model (RFIM). As shown by 
Cardy and Jacobsen& in the limit of d — > 2 and g — *■ cxd 
the interface Hamiltonian of the two problems have the 
same type of solid-on-solid (SOS) description, from which 
follows that the energy exponents of the tricritical RBPM 



are equivalent to the magnetization exponents of the crit- 
ical RFIM. Furthermore, analyzing the rcnormalization 
group (RG) flow it was conjectured that the above map- 
ping stays valid for d > 2, in particular at d = 3. Thus 
the tricritical exponents should be q independent, at least 
for large q, leading to the same exponents for any disorder 
induced tricritical points. 

These conjectures, which could be of experimental rel- 
evance, have not yet been verified numerically. The in- 
accuracies in the simulations have mainly two sources: 
i) it is difficult to precisely locate the tricritical point 
due to strong cross-over effects and ii) for not too large 
q the tricritical disorder is quite small which results in 
large breaking-up lengths in the systeniiS, thus one has 
to treat quite large lattices. 

In the present paper we consider the 3d ferromagnetic 
Potts modeli^ for large value of q, which has a strongly 
first-order transition for non-random couplings and study 
the effect of bond randomness including also the case of 
bond dilution. In the large-g limit the high-tcmpcrature 
expansion of this problem^^ is dominated by a single 
diagramii which is exactly calculated by a combinatorial 
optimization methodic. Weak disorder is studied pertur- 
batively, whereas for stronger disorder we have performed 
extensive numerical calculations. In particular we have 
studied the properties of the tricritical point and checked 
its possible relation with the critical fixed point of the 
RFIM. Wc have also studied the form and universality, 
i.e. disorder independence, of the critical singularities. 
For theses quantities some results have already been an- 
nounced in a letteria. 

The structure of the paper is the following. The RBPM 
and the optimization method used in the study for large 
q is presented in Scc^ Perturbative treatment of the 
problem in the weak disorder limit is shown in Sec lIIII 
Numerical results about the phase diagram as well as 



2 



on critical and tricritical singularities are given in Sec II VI 
and discussed in SeclVl 

II. MODEL AND ITS PHASE DIAGRAM 

The g-state Potts modeli^ is defined by the Hamilto- 
nian: 

H = -^JijS{ai,aj) (1) 

in terms of the Potts-spin variables, ai = 0, 1, • • • , g — 1. 
Here i and j are sites of a cubic lattice and the summa- 
tion runs over nearest neighbors. The couplings, > 0, 
are ferromagnetic and identically and independently dis- 
tributed random variables. In this paper we use a bi- 
modal distribution: 

P{J,,)=pSiJil + S)^J,,) + {l-p)S{Jil-S)-J,j) (2) 

with p = 1/2. The parameter < 5 < 1 plays the role of 
the strength of disorder and ai S = and ^ = 1 we recover 
the non-random and the diluted systems, respectively. 

For a given set of couplings the partition function of 
the system is conveniently to write in the random cluster 
representationi^ as: 

z = Y.i''''^ n b''--^ (3) 

G ijeG 

where the sum runs over all subset of bonds, G and c{G) 
stands for the number of connected components of G. In 
Eq.Q we use the reduced temperature, T ^ Tlnq and 
its inverse (3 — > f3/hiq, which are of 0(1) even in the 
large-g limiljiS,. In this limit we have g^-^'j ^ 1 and the 
partition function can be written as 

Z=Y. I^^""^ ' '^(G) = c(G) +(3Y. (4) 

GCB ijGG 

which is dominated by the largest term, (p* = maxc (t>iG). 
Note that the optimal set itself generally depends on the 
temperature. The free-energy per site is proportional to 
(j)* and given by — /?/ = (p* /N where N stands for the 
number of sites of the lattice. 

The optimization problem in Eq.l^J contains a cost- 
function, 0(G), which is sub-modulai— and there is an 
efficient combinatorial optimization algorithm, which at 
any temperature, works in strongly polynomial timai^. 
This algorithm finds a set of bonds which minimizes 
the cost- function. We call such a set an optimal set. 
The variation of the optimal set with the temperature 
is illustratedSi in FiglJ for S = 0.875 and L = 24. At 
low temperature the optimal graph is compact and the 
largest connected subgraph contains a finite fraction of 
the sites. In the other limit, for high temperature, most 
of the sites in the optimal set are isolated and the con- 
nected clusters have a finite extent, the typical size of 
which is used to define the correlation length, ^. 



Between the low-tempcrature (ordered) and the high- 
temperature (disordered) phases in the thermodynamic 
limit there is a sharp phase transition. The numerically 
calculated phase-diagram as a function of the temper- 
ature, T, and the disorder, 5, is shown in Fig. |3 A 
detailed analysis of the phase-diagram is postponed to 
Sec lIVI Here we just note that the transition is of first 
order for weak disorder, in which case ^ stays finite at the 
transition point, but the transition is of second order for 
strong enough disorder, when the correlation length is di- 
vergent at the transition point. This second possibility is 
illustrated in the central part of Fig^ when at the tran- 
sition point the largest connected cluster is a fractal and 
its fractal dimension, df, is related to the magnetization 
critical exponent, see in Sec lIVI The fractal dimension 
of the giant connected cluster at the transition point is 
shown in the inset of Fig. [3 as calculated in Sec lIV Cl 
In the first-order regime it is df = 3, i.e. the cluster is 
compact, whereas in the second-order regime it is d/ < 3 
and practically independent of the strength of disorder. 
At the tricritical disorder, St, the giant cluster is a frac- 
tal but its fractal dimension is different from that in the 
critical regime, see in Scc lIVDl 



III. PERTURBATIVE CALCULATION FOR 
WEAK DISORDER 

A. Non-random model and Imry-Ma type 
argument for weak disorder 

In the non-random model, 6 = 0, there are only two 
homogeneous optimal sets, which correspond to the T = 
and T ^ oo solutions, respectively, see in RefiSi. For 
T < Tc(0) it is the fully connected diagram with a free- 
energy: -PNf = 14- Npjd and for T > TdO) it is 
the empty diagram with —/SNf = N. Consequently the 
transition point is located at: rc(0) = Jd/{1 — 1/N) and 
the latent heat is Ae/Tc(0) = 1 - 1/A^. 

In the presence of disorder, 6 > 0, the optimal sets are 
homogeneous only in a limited temperature range and 
new non-homogeneous optimal diagrams appear, see in 
Fig^ In the disordered phase in the vicinity of Tc{0) 
the typical linear size of the connected clusters, /, can 
be estimated along the lines of the Imry-Ma argumenliS^. 
Adding to the optimal set a cluster with a number of 
sites of the order l'^ will decrease the number of connected 
components, and therefore increase the cost-function by 
l"^. On the other hand, it will add a number of the order 
dl"^ — dl'^'^^ of bonds, each having a weight (1 ± S)/d 
since the temperature is close to Tc{0). It results in a 
competition between a term behaving like Z''/^, which 
represents the gain due to disorder fluctuations and a 
term like Z''"^, which is the lost due to the creation of an 
interface. In d = 2 the two terms balance each other and 
extreme fluctuations of the disorder will create clusters 
of unlimited sizeifl. As a consequence in 2d there is no 
coexistence of pure phases and the transition is of second 
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T= 3.200 



FIG. 1: Connected parts of typical optimal sets with L = 24 
at 5 = 0.875. Top: T = 3.122J < T^, centre: T = 3.141J ^ Tc 
and bottom: T = 3.200J > Tc. 
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FIG. 2: Phase diagram of the random-bond Potts model 
for bimodal disorder. In the ordered phase in the optimal 
set the non-connected sites are percolating for 5 > Spr and 
Tpr{5) < T < Tc{5). Numerical results indicate, that the 
first- and second-order transition regimes are separated by 
a tricritical disorder, 5c, which corresponds to the border of 
the percolation regime, Sc — Spr. Inset: fractal dimension of 
the giant connected cluster at the transition point as calcu- 
lated in Sections IIV CI and IIV PI The straight lines at df — 3 
and df = 2.40 indicating the values at the first-order and the 
second-order transition regimes, respectively, are guide to the 
eye. 



order for any small amount of disorder—. 

On the contrary in d > 2, in particular in three dimen- 
sions for weak disorder the surface term is dominating, 
thus the connected clusters have a finite extent and the 
transition is of first-order in accordance with the phase- 
diagram in Fig. 12 In this case connected clusters are 
created due to extreme fluctuations of disorder, and the 
statistics of these rare regions will be considered in the 
following subsection. A similar analysis of the ordered 
phase will be given afterward in Sec. 1111 CI 



B. Disordered phase 

Here we make the estimate of the free-energy as simple 
as possible and therefore we consider rare regions of the 
shape of cubes of linear linear size, Z > 2, in which all 
the n+ = 3(/'^ — P) internal coupHngs are strong, being 
J{1 + S). The cost-function of the diagram in which all 
the bonds of the cube are present relative to the empty 
graph is given by: 



Af+{1) = 3/3(1^ -f)J{l + S)-{l^-l) 



(5) 



The optimal set is then the inhomogeneous one contain- 
ing the connected cluster, if A/+(/) > 0, which could 
take place if the disorder satisfies the relation, 5 > S^{1), 
with 



(6) 
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In another words for a given disorder, 5, there is a hmiting 
size, l+{5) = (1 + Vl + 4:6)/25 « 1/S, and only large 
enough clusters with I > l+{5) can exist in the optimal 
graph, provided the temperature is in the range of 3 < 
T/J < T+{l)/J = 3 + 3((5 - + 5+{l)). Thus for 

weak disorder only large (and very rare) clusters can be 
found and as the disorder is increased at discrete values 
of 5 new, smaller and more probable connected clusters 
will appear. This mechanism will lead to discontinuous 
behavior of the free energy as a function of the disorder, 
which is indeed observed in numerical calculations. 

The probability of having all the couplings strong 
in a cube is exponentially small, P+{1) — 2~"+, thus 
the density of connected cubes of size I in the optimal 
set is given by p+{l) = 2~'^'' ~' \ The free energy of 
the non-homogeneous optimal set is given by the sum of 
contributions of the connected clusters of different size: 

-PNU^N + N J2 A/+(Op+(0. (7) 

1>1+{S) 

Since p+{l) is a very rapidly decreasing function of / the 
sum in Eq.jTI) is dominated by the term with I = 1+{S). 



C. Ordered phase 

In the ordered phase, T < Tc(0), we consider rare re- 
gions of extreme disorder fluctuations also of the shapes 
of cubes with / > 1, so that all the couplings, n_(Z) = 
+ P), starting from the points of the cubes are weak, 
being J{l — S). The cost-function of the diagram in which 
all weak bonds of the cube are absent, i.e. there is a cube 
of isolated points embedded into the full diagram, rela- 
tive to the full graph is given by: 

A/_ (0 = /3 „ 3^(^3 _ ^) ^g) 

As for the disordered phase the optimal set is the inhomo- 
geneous one containing the isolated points, if A/_ (1) > 0, 
which is the case for strong enough disorder, 5 > 6-{l), 
with 

MO-^. (9) 

Now for a given disorder, 6, the limiting size is, = 
1/(5 — 1 « 1/S, and the I > l-{6) "empty clusters" exist 
in the optimal graph in the range of temperature 3 > 
T/J > T^{l)/J = 3 -h 3(5 - 5-{l))/{l + S-{1)). Now the 
density of cubes of isolated points of size I in the optimal 
set is given by p-{l) = 2"^^' ' and the free energy of 
the non-homogeneous optimal set can be written as: 

- ^Nf_c^l + Npjd + N A/_(Op-(0, (10) 

1>1-{S) 

which is dominated by the term with I — 1-{S). 



D. Phase transition 

The phase-transition temperature, Tc{d), is obtained 
by equating the free-energy in the two phases: /_ = /+. 
Since it is comparatively easier to create an "empty clus- 
ter" in the ordered phase, than a connected cluster in the 
disordered phase which results in the shift of the phase 
boundary towards higher temperatures as the disorder 
is increasing, see Fig. ^ Further observation, that the 
phase-boundary is discontinuous at 5+(/) and S-{1), for 
all integer l-s. For small (5, which corresponds to large 
l-s these jumps are very frequent and in a mathemati- 
cal point of view the phase boundary is a non-analytical 
function of S. However these jumps are very small and 
the phase-boundary can be approximated by a continu- 
ous curve which asymptotically behaves as: 

ln(T,-T,(0))~-l. (11) 

We obtain similarly for the latent heat and for the jump 
of the magnetization at the transition point: 

\n{Ae/Tc-l) ^\n{Am-l) ^ . (12) 

Thus the phase transition stays first order and as the 
disorder is switched on there is an essential singularity in 
the thermodynamical quantities as a function of S'^. 

E. Distribution of the finite-size transition 
temperature 

For a large finite system of linear size, L, and for a 
given realization of the disorder one can define a transi- 
tion temperature, Tc{L), at which the disordered phase 
and the ordered phase (having a spanning giant cluster) 
coexists. The distribution of Tc(L) is discontinuous for a 
finite system but expected to be Gaussian for large L, so 
that characterized by its average, T^'"{L), and its vari- 
ance, var[Tc{L)] = [ATc(L)]^. The shift of the average is 
asymptotically given by: 

^r(i)"^c(oo)~L-l/^ (13) 

whereas the width scales with another exponent, as: 

ATc{L) - L-'^/" . (14) 

These relations define the exponents and and hold 
for second orderi^SiSiiSSiSSiSi as well as for first order 
phase transitionsS^iSS. At a first-order transition what 
we study here i) is just the discontinuity fixed point 
value of the correlation length exponenli^fliSi and given 
hy D = 1/(1 = l/d. This exponent describes the varia- 
tion of a diverging length-scale which can be measured 
in a typical sample. The width of the distribution can be 
obtained from the consideration that the density of clus- 
ters, p±, has a fluctuation in finite systems which has a 
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TABLE I: Breaking-up disorder in the ordered phase: com- 
parison of the interpolation formula in Eq. 1151 with numerical 
results 



L 


"int 


5^ 


7 


0.466 


0.469 
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0.454 


0.458 


15 


0.434 


0.438 


17 


0.430 


0.435 


23 


0.420 


0.429 


31 


0.412 


0.423 


39 


0.406 


0.418 



width of, Ap±{L) / p± ^ pL^'^^^, according to the central 
limit theorem. The fluctuations in Ap±{L) are then seen 
in the fluctuations of Tc{L), too, leading to an exponent, 
1/ = 2/d. This exponent is related to another diverging 
length-scale, which can be deduced from the study of av- 
erage quantities. Note that this length is not present in 
the pure system and by switching on disorder its pref- 
actor, p±, presents an essential singularity as a function 
of 6. These results are in accordance with the considera- 
tions presented in Sec. VII of Ref . ■ Numerical analysis 
of Tc{L) in for stronger disorder is given in Sec lIVI 



F. Breaking-up lengths 

For weak disorder the density of elementary excitations 
is very small and the average distance between two exci- 
tations is given by ~ and L5_ ~ p-{S)~^/^, 
in the two phases, respectively. and L'i_ can be in- 
terpreted as breaking-up lengths, since in a finite system 
of linear size, L < Lj_, the optimal set is homogeneous. 
Using results about the critical densities we obtain: 

log2 L\ = {l + 2S- 25^ + VTTl5)/25^, 

log., Lt = {S-'- if + {S-'-l)\ (15) 

which in principle is valid only at i5 = 6±{l), however 
as an interpolation formula we can use them for not too 
small ^-s, too. 

It is easy to see that the braking-up lengths in the or- 
dered phase are much smaller, then in the disordered one. 
For example at the value of 6 = 1/2, which corresponds 
to the existence of a single hole in the ordered phase and 
to a ^ = 2 connected cube in the disordered phase, re- 
spectively, L'L = 4, whereas i'j. ~ 2^^^. We have checked 
the accuracy of the formula in Eq. H15() by comparing the 
predicted breaking-up disorder, S'L, for a given size, L, 
with that calculated numerically. As seen in TablcUthere 
is a satisfactory agreement, even for not too small S-a. 



IV. NUMERICAL RESULTS 

In the numerical calculation we have treated sam- 
ples with random couplings having a cubic shape with 
periodic boundary conditions and a linear size L = 
16, 24, 32. For S = 0.875 we went up to L = 40 and 
in some cases we made calculations for odd values of 
L, too. The free-energy as well as the magnetization 
is calculated exactly by the combinatorial optimization 
algorithm, called as "optimal cooperation"— and aver- 
aging is performed over several thousands of samples, for 
the largest size the number of realizations was several 
hundreds. For a fixed temperature the optimal coopera- 
tion algorithm works in strongly polynomial time. As we 
have already mentioned in the application of the method 
in 2(iiS for a finite system the free-energy is a piece-wise 
linear function of the temperature. For the bimodal dis- 
order the number of linear parts, Np, is found to increase 
with the size of the system L as well as with the disorder 
strength 6. A rough estimation yields 

Np ~ CL^+^ , (16) 

with the constant C around one half. We have managed 
to implement our method in such a way that we can 
calculate the free-energy in the whole temperature range, 
i.e. in all linear parts. Note that with the hypothesis 
((TE|l the exact calculation of the free energy in the whole 
temperature range is still polynomial, whatever the value 
of the disorder S, 

As we have shown in Sec lIIIFl the breaking-up length 
is large for weak disorder and it is of the order of L'' = 
40 for i5 w 0.4. Therefore we have restricted ourselves 
to the disorder range: 0.4 < S < 1, which contains all 
interesting parts of the phase diagram. 



A. Phase diagram 

The numerically calculated phase diagram as a func- 
tion of disorder and temperature is already presented 
in Scc|n] in Fig. |21 here we give a detailed analysis of 
the results. The phase boundary between the ordered 
and disordered phases is calculated by considering two 
quantities. Finite-size or percolation transition temper- 
atures, Tc{L), are identified in each sample at the point 
where the largest connected cluster of the optimal set 
starts to percolate the finite sample. The distribution 
of Tc{L) is shown in Fig|31 in the first-order transition 
regime {S = 0.5) and in Fig^in the second-order transi- 
tion regime {S ~ 0.875). 

While for d = 0.875 the distribution of the transi- 
tion temperatures is well described by a Gaussian even 
for relatively small systems, for S ~ 0.5 the distribu- 
tions show deviations from a Gaussian. As we can see in 
FigElin the first-order regime the distribution for finite 
L consists of well separated peaks and it becomes (quasi- 
)continuous only in the thermodynamic limit. Also the 
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FIG. 3: Distribution of the finite-size percolation transition 
temperatures in the first-order transition regime at S — .5 
for three different sizes. The Gaussian distributions having 
the same average and variance are normalized and serve as a 
guide to the eye.. 



distributions are non-symmetric and the skewness seems 
to vanish slowly with the size of the system. It is also 
evident from FigsOlandElthat the finite-size shift of the 
average transition temperature has a different sign in the 
two regimes and it is much smaller for a first-order transi- 
tion in particular if we compare with the variation of the 
width of the distributions. On the contrary in the second- 
order transition regime these two characteristics of the 
distribution arc in the same order, see in Fig^J These 
observations are in accordance with the scaling picture 
in Sec. IIII El which predicts two different exponents, 
and V in the first-order transition regime, whereas in the 
second-order regime these exponents are the same. How- 
ever, due to the comparatively small sizes we have and 
the deviations of the distributions from the Gaussian we 
could not numerically estimate these exponents. 

We have also calculated the transition temperature 
from the position of the maximum of the averaged specific 
heat, Cy{T,L), which is expected to be shifted in finite 
systems as the maxima of Tc{L), as given in Eq.lQTJ. In- 
deed the numerical results on the specific heat in Fig. 
at (5 = 0.875 are compatible with the same estimate for 
Tc, as obtained in Fig^ through the percolation transi- 
tion temperatures. 

As discussed before magnetic order in the system is 
related to the properties of the largest connected cluster 
in the optimal set: it has a finite extent in the disordered 
phase, whereas a finite fraction of sites belongs to this 
cluster in the ordered phase. The optimal set in the or- 
dered phase, however, can be of two different kinds as 
far as the structure of isolated points is considered. For 
weak disorder, i5 < Spr, or for low temperature, T < Tp^, 
the isolated sites form finite clusters. This is always the 
case in 2d. In 3d, however, and this is a new feature of 
these systems for strong enough disorder and high enough 



FIG. 4: The same as in Fig^jin the second-order transition 
regime at S — 0.875 



temperatures the isolated sites percolate the sample, too. 
This percolation regime is also indicated in Fig|21and we 
argue below that its existence close to Tc is necessary 
to have a second-order transition in the system. Indeed, 
the correlation length in the ordered phase is given by 
the size of the largest connected finite cluster, which is 
isolated from the giant connected cluster. Since this large 
finite cluster is embedded into isolated points its size can 
be divergent at Tc only if the isolated sites percolate the 
sample. In this way we obtain for the value of the tricrit- 
ical disorder, St, separating the first- and second-order 
transition regimes that it satisfies the relation, 

6t > Spr . (17) 

Numerical results which arc presented below are in favor 
of the conjecture that in this relation the equality holds. 

For the percolating transition temperature of the iso- 
lated sites, Tpr{S), we make the following calculation. In 
the dilute model, S = 1, all the strong bonds, Ji = 2 J, 
are present in the optimal set, provided the temperature 
is below the value of Ji. For T > Ji, however, in the op- 
timal set there are no dangling bonds, i.e. by increasing 
the temperature over T = Ji sites which have just one 
strong bond are removed from the optimal set. Since the 
dangling bonds are a finite fraction of the bonds^ the 
non-connected sites become percolating at Tpr(l) = 2 J. 
For S < 1 there are also weak J2 bonds in the system 
and the removal of one strong dangling bond from the 
optimal set is accompanied by the removal of some weak 
bonds at the same time. The average number of removed 
decorating weak bonds is four which is possible in the 
temperature range: 

Tc>T> Tpr{5) = ,/i + 4J2 = ./(5 - 3(5) . (18) 

The numerical results show that at Tpr{S) in a finite frac- 
tion of samples there is a giant cluster of isolated points 
which spans the finite cube, see in Fig|31 
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FIG. 5: Magnetization and percolation probabilities of the 
optimal set ("Present": for the occupied bonds of the giant 
cluster; "Absent" : for the isolated sites) of the random model 
with S = 0.875, as a function of the temperature in a finite 
lattice with L = 24. The arrow at T = Tpr = 2.375J indi- 
cates the percolation temperature in Eg. 11181 and for T < Tpr 
the magnetization is close to that for ordinary percolation. 
The singular jumps in the magnetization are a consequence 
of the discrete form of the probability distribution. Inset: in 
the vicinity of the transition point the magnetization and the 
percolation probability of the present bonds are close to each 
other, the differences are due to finite-size effects. 

At the percolation transition temperature there is a 
sudden change in the structure of the giant connected 
cluster, which results in singularities in the thcrmody- 
namical quantities. As an illustration we show in Fig. 
|S1 the temperature dependence of the magnetization at 
S = 0.875, i.e. in the second-order transition regime. The 
magnetization goes to zero at the phase-transition point, 
which - in the thermodynamic limit - coincides with the 
percolation transition of the occupied bonds. For finite 
systems the magnetization and the probability of having 
a spanning cluster has small deviations, as illustrated in 
the inset of Fig. |S1 The magnetization has another sin- 
gular jumps for T < Tc, which are due to the discrete 
form of the disorder. Among these singularities the most 
pronounced is that around the percolation temperature 
at Tpr = 19/8. As seen in Fig. Othc percolation prob- 
ability of the isolated sites has a finite value at Tpr and 
it goes to the value of one within a small temperature 
range. It is expected that by decreasing the disorder to 
the tricritical value, the singularities at Tc and Tpr merge 
into a new type of tricritical singularity. 



B. Order of the transition 

To decide about the order of the transition one gener- 
ally studies the behavior of the latent heat in the system. 
This type of analysis, however, for the large-g state Potts 
model is complicated, if the disorder is discrete, as in our 
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FIG. 6: Temperature dependence of the magnetization in the 
transition region for L = 16 with different strength of disor- 
der. From left to right 5 = 0.5 - first-order regime; St = 0.6593 
- tricritical point; 5 = 0.925 - second-order regime. 

case. As already discussed in 2d~^ the internal energy 
of the system displays discontinuities, both at and out- 
side the critical point. These discontinuities, however, 
are generally connected to such degeneracies of the opti- 
mal set which are related to the removal of finite num- 
ber of bonds, thus do not modify the global structure of 
the giant connected cluster. As a consequence these sin- 
gularities can not be interpreted as a sign of first-order 
transition. In 3d and for the bimodal disorder used in 
this paper this type of non-generic discontinuities of the 
internal energy are also present therefore we did not try 
to make an analysis of this quantity. 

Instead we have studied the magnetization in the sys- 
tem, which is defined as the fraction of sites in the largest 
connected cluster. A jump in the magnetization indi- 
cates a fundamental change in the shape of the largest 
cluster, thus a first-order transition. The magnetization 
as a function of temperature is shown in Fig [3 for differ- 
ent values of the disorder corresponding to the different 
transition regimes. A direct measurement of the jump in 
the magnetization is possible only up to (5 < 0.6, see in 
Ref^, whereas the lower limit of the second-order tran- 
sition regime is estimated through the calculation of the 
fractal dimension of the giant cluster at the transition 
point, see in Fig|21 This type of analysis at Spr (see in 
Sec. IIVD|I has given a fractal dimension, d*^ < 3, which 
together with Ea. H17|) indicates that the first-order tran- 
sition regime stops at Spr- 



C. Second-order transition regime 

In the second-order transition regime we have 
made detailed calculations at the disorder 6 = 
0.75, 0.8, 0.875, 0.925 and at J = 1. The most de- 
tailed studies are performed a,t 6 ~ 0.875, in which case 
the largest system is L = 40 for the other cases we went 
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up to L = 32. With these investigations our aim was to 
check universahty of the critical properties of the system. 

The magnetization exponent, /3, and the magnetiza- 
tion scaling exponent, x = P/v, is related to the frac- 
tal dimension of the giant cluster, df, as d — a; = df. 
Among these critical parameters it is the fractal dimen- 
sion which can be determined with the highest preci- 
sion. To obtain the fractal dimension we considered a 
reference point of the percolating cluster and measured 
the mass (number of points) in a shell around the refer- 
ence point with unit width and radius, r. The average 
mass, s{r,L,T), is expected to scale close to the transi- 
tion point: t = {T - Tc)/Tc < 1 as: 



sir,L,t) = L'^f-^§{r/L,tL^/'') 



(19) 



Generally, one sets the second argument of the scaling 
function, s{p,t), to be zero by performing the calcula- 
tion at the transition temperature, T = Tc- In our case 
Tc is not known exactly, therefore we have used another 
strategy. For each size we set, T = Tc{L), which is the 
average finite-size (percolation) temperature at the given 
size. With this choice the second argument of the scaling 
function, t is asymptotically constant and thus the scal- 
ing function depends only on one parameter: s = s{r / L). 
Our scaling picture is checked in Fig. [7| in which the scal- 
ing plot of the mass in the shell is shown for 5 = 0.875. 
The accuracy of the scaling collapse is measured by the 
area of the collapse region which is shown in the inset of 
Fig. 13 It is seen that the optimal scaling collapse is ob- 
tained with a fractal dimension, dj — 2.40(2). This type 
of analysis of the fractal dimension is repeated for an- 
other values of the disorder, too. Since in these cases the 
available sizes of the systems are comparatively smaller 
we have somewhat larger errors. The fractal dimensions 
are shown in the inset of Fig. |21 

The correlation length critical exponent, v, can be cal- 
culated from the shift of the finite-size critical temper- 
ature, as given in Ea. ljlHII . However this method has 
quite a large error. A more accurate estimate can be ob- 
tained from the scaling behavior of the magnetization: 
m(t, L) = L^'"rh{tL-^/'^). Here we set Xm = d — df from 
the previous calculation and from an optimal scaling col- 
lapse as shown in Fig. |Slwe have obtained v = .73(2). 
(The area of the collapse region as a function of i/ is 
shown in the inset of Fig. |H|) Note that v satisfies the 
rigorous bound for disordered systems^: v > 2/d. 

Finally, we have investigated the behavior of the spe- 
cific heat at the transition point. As seen in Fig. |51the 
maximum of the specific heat is increasing with the size, 
therefore from the finite-size scaling result: C^'"^(L) ~ 
X"/" one would conclude a > 0. This is, however, in 
conflict with hyperscaling and with the bound of v in 
disordered systems. Therefore we tried to fit the numer- 
ical date by including a constant in the r.h.s. of the 
finite-size scaling form. The fit in this way, however, is 
not satisfactory. In order to get a non-positive a, one 
should have a constant, which is more than one order of 
magnitude larger, than the finite-size data. Therefore we 




FIG. 7: Scaling plot of the mass of a shell of the infinite 
cluster, for each size at the average spanning temperature, 
see text, at S — 0.875. Optimal collapse is obtained with 
a fractal dimension, df = 2.40. Inset: area of the collapse 
region as a function of df. 




FIG. 8: Scaling plot of the magnetization at 5 = 0.875 as 
a function of the distance from the critical temperature. By 
fixing P /v = d — df = 2.4 the best collapse is obtained with 
V — .73(2). In the inset the area of the collapse region is 
shown as a function of u, the arrow indicates the position of 
the minima. 



have concluded that the asymptotic regime of the spe- 
cific heat is very far from the possibilities of present day 
numerical calculations. 



D. Tricritical transition 

Analyzing the structure of the optimal set in the or- 
dered phase we have got a relation between the tricrit- 
ical disorder and the percolating disorder, as written in 
Ea. (|17|) . Later we have conjectured that in this relation 
probably the equality holds, i.e. at the tricritical point 
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FIG. 9: Finite-size specific heat as a function of temperature 
at (5 = 0.875. According to liyperscaling the maxima of the 
curves should approach a finite limiting value. The relatively 
large errors are due to numerical derivation of a piece-wise 
linear function 



three regions (disordered phase, ordered phase and per- 
colating regime) meet. Here we check this conjecture nu- 
merically in the following way. We sit on the analytical 
continuation of the percolation transition line in Eq. I|18() 
for 6 < Spr (thus for T > Tpr) and calculate the finite- 
size transition point, which is located at such 5t{L) for 
which in half of the samples the giant connected clus- 
ter percolates the finite cube. Then we have analyzed 
the fractal properties of the giant connected cluster and 
repeated the procedure as described in Sec lIVCI for the 
second-order transition regime. Using the relation about 
the number of sites in a shell in Ea. (|19|l and fixing the 
temperature for each size through St{L) we have made a 
scaling analysis as in FigElfor the second-order transition 
regime. According to the results in Fig^Jthe giant clus- 
ter is a fractal and the best collapse of data is obtained 
(see inset) with = 2.90(2). This value is definitely dif- 
ferent from that at a first-order transition, df = d = 3, 
and also diifers from that in the second-order transition 
regime. Thus our numerical results are in accordance 
with the conjecture, that the tricritical disorder is given 

by 



Our numerical data are very sensitive to the value of 
(5, and by extrapolating the finite-size results obtained 
from the best collapse we find S = 0.65930(5) and 
T = 3.0221(1). Furthermore, we obtain the estimate for 
the anomalous dimension of the tricritical magnetization, 
0.10(2). 

In order to calculate the correlation length exponent, 
v, at the tricritical point we consider the magnetiza- 
tion, m{t,T, L), as a function of reduced temperature, 
t, difference of the tricritical disorder, t = [5t — S)/5t, 
and the size, L. According to finite-size scaling we 



FIG. 10: Scaling plot of the mass of a shell of the infinite 
cluster analyzed along the percolation transition temperature, 
Eq.lESJ. For each size we fix the disorder and thus the tem- 
perature as the average spanning temperature, see text. Op- 
timal collapse is obtained with a fractal dimension, dj^ — 2.94. 
Inset; area of the collapse region as a function of dj. 



have m{t,T,L) = L^^^^mirL^/^^ ,tL^/^) and by having 
a disorder dt{L), (thus t{L)) for each finite size we fix 
asymptotically the first argument of the scaling func- 
tion. Consequently form the optimal scaling collapse of 
m{t, t{L), L)L^'-^ we can obtain the correlation length ex- 
ponent. We noticed that at the tricritical point the esti- 
mate for u is more sensitive to the range in which the col- 
lapse is performed. Using the symmetric region indicated 
in Fig. 111! the minimum of the scaling area in the inset 
of Fig. ^]is aX V = 0.64. However using an asymmetric 
window, -0.1 < (T - Tc)Li/'' < 0.3, we obtain a value 
which is slightly larger than the borderline value u = 2/2i 
according to the criterion by Chayes et alS. Therefore we 
conclude the estimate of the correlation length exponent 
at the tricritical point as: v = 0.67(4). 

The tricritical fixed point of the random bond Potts 
model in the large-g limit is related to the critical fixed 
point of the random-field Ising model according to a 
mapping due to Cardy and Jacobsen^. As described in 
Ref.jg the interface Hamiltonian of the two problems in 
the solid-on-solid approximation are equivalent to each 
other in d = 2 + e dimension, and this mapping is ex- 
pected to hold for larger values of d, in particular for 
d = 3. According to the mapping magnetization-like ex- 
citations of the RFIM correspond to energy-like excita- 
tions in the random bond Potts model. In particular the 
correlation-length exponent at the tricritical point of the 
RBPM is conjectured to relate on the critical exponents 
of the RFIM as 



i. = i.«^/(/?«^ + 7«^) . (20) 



For the 3d RFIM with Gaussian disorder recent estimates 



10 




-0.4 -0.3 -0.2 



(T-TJ L" 



FIG. 11: Scaling plot of the magnetization at the finite-size 
tricritical disorder as a function of the distance from the tri- 
critical temperature, see text. By fixing (i^ /u = d — d^f — 2.90 
the minimum of the collapse region (see inset) is obtained 
with V — .64. 



on the critical exponents ar( 



,RF 



1.37(9), /3^^ = 0.017(5), 7 



RF 



2.04(14) 



(21) 

which imply through Ea. H20() ~ 0.67(2), which coin- 
cides with our direct calculation. In this way our study of 
the tricritical singularities lends support to the mapping 
by Cardy and Jacobsen. The tricritical magnetization 
exponent, a;^j, we calculated above can not be predicted 
by the mapping therefore it gives a completely new piece 
of information. 

Finally, we note that in the spirit of the RG de- 
scription the number of Potts-states, more precisely 
In (7, is a dangerous irrelevant scaling variable. There- 
fore the calculated tricritical exponents are expected to 
be g-independent and probably universal for any three- 
dimensional systems in which the first-order transition 
in the pure systems is soften by tricritical disorder. As a 
matter of fact, it is much easier to compute the free en- 
ergy of the Potts model at q infinite, using submodular 
functions theory, than at q very large. This will make dif- 
ficult to check the assumption above, at least numerically 
and in the framework of the Potts model. 



outside the critical point and this regime is called the 
disordered and ordered Griffiths phasea^^. The singular 
behavior is due to rare regions which, due to strong disor- 
der fluctuations, are locally in the non-stable thermody- 
namic phase of the system. It is easy to see that similar 
effects can be observed outside and at a first-order transi- 
tion point, as we show in the following for our model with 
weak disorder. Note, however, that our numerical results 
concern q infinite, for which no dynamics is defined. In- 
deed, only the spin-state compatible with the diagram(s) 
maximizing the function in Eq.Q are possible, all these 
states being equiprobable. The analysis below refers to 
q very large, but finite. 

In the disordered Griffiths phase of the random bond 
Potts model such a rare region is a domain (cube) of lin- 
ear size I, having only strong couplings, as described in 
Sec lIII B1 These regions are indeed rare since they ap- 
pear with a probability: P{1) ^ 2'^(' ~^ \ In such a clus- 
ter in equilibrium all the spins are typically in the same 
state, provided the conditions in Sec lIII Bl are fulfilled. 
However during a relaxation process they flip into an- 
other parallel configuration. Having heat-bath dynamics 
the relaxation time, tr, can be estimated in such a way, 
that in a thermally activated process the cluster should 
overcome an energy barrier, AE(l), which is the energy 
of creation an interface in the system. The relaxation 
time is then given by: tr ^ exp{/3AE{l)). Generally the 
interfacial energy is proportional to the number of sites 
involved in the interface: A{1) ~ l'^~^£,±, where is the 
width of the interface. It is knowni2& that outside the 
first-order transition point is finite, whereas at the 
transition point it is divergent as ^± ~ with a ther- 
mal wandering exponent: C = (3 — d)/2. In particular 
at d = 3 we have a logarithmic divergence. Now cal- 
culating the distribution of relaxation times we obtain 
P{tr) ~ exp[— est [IntrY/^), where a; = — 1, outside 
the transition point of the cluster and lo = d — 1 + Ci close 
to the transition point. Then the average autocorrelation 
function is given by: 



G(t) 



dtrP{tr) exp{—t / tr) ^ cxp —cst (Int) 



\d/u 



.(22) 

which has a different form at the first-order transition 
point and in the disordered phase. 



DISCUSSION 



Effect of the form of disorder 



Here we discuss some aspects and possible extensions 
of our results obtained in the previous Sections. 



A. Dynamical behavior for weak disorder 

In a random ferromagnetic system having a second- 
order transition point different dynamical quantities (sus- 
ceptibility, autocorrelation function, etc.) arc singular 



In this paper in the numerical studies we considered 
bimodal i.e. discrete form of disorder. This was mainly 
due to technical simplifications. The basic behavior of 
the system in particular its critical properties are not ex- 
pected to change using a continuous form of disorder. 
For example the phase-diagram has the same topology 
as shown in Fig|21 i.e. disordered and ordered phases 
and first- and second-order transition regions. The per- 
colation region of the ordered phase also exists and very 
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probably the tricritical point is located at the meeting 
point of these three phases or regions. Also the critical 
singularities at the second-order transition line are very 
probably disorder independent, as we have already ob- 
served in the two-dimensional problemifl. To the same 
question for the tricritical transition is more difficult to 
answer. Having in mind that the 3d RFIM might have 
disorder dependent singularities^ the same can be true 
for the tricritical singularities, due to the mapping as de- 
scribed in Sec II V Ul Finally the discrete nature of the 
disorder can result in non-physical discontinuities in the 
internal energy, which are washed out by continuous dis- 
order. 



C. Effect of the number of states - q 

In the second-order transition regime the critical ex- 
ponents are g-dependent, which can be seen from the 
results of numerical investigationsi^iii and the same sce- 
nario holds in two dimensions, too. The tricritical sin- 
gularities, however, at least those which are related to 
the energy density, are very probably g-indcpendent and 
thus 'hyperuniversaP. It would be very interesting to 
check this statement for another systems, since the Potts 
model in three dimensions for q finite but very large, 
seems out of the reach of the usual methods. Finally 
we mention that the non-physical discontinuities in the 
internal energy are also absent for finite value of q. 

D. Model with correlated disorder 

Finally, we consider our model with correlated disor- 
der, in which translational invariance is present in the 
vertical direction, in which the couplings are constant, 
J_L, whereas in the horizontal 2d planes the couplings 
arc random, Jy , and strictly correlated in each planes. 



This type of columnar disorder is introduced by McCoy 
and Wu in the 2d Ising model^. Due to translational 
symmetry in the vertical direction the model is conve- 
niently studied in the transfer matrix formalism. Us- 
ing the extreme anisotropic limit of the model^, when, 
J_\_/ Jij — > the transfer matrix is written into the form: 
T = exp(— rTi), where r is the infinitesimal lattice spac- 
ing and Ti is the Hamiltonian operator of the 2d quan- 
tum Potts model with random couplings. This latter 
model can be studied by a strong disorder renormal- 
ization group method^ which leads to g-independcnt 
critical properties. According to numerical results^Siii 
the correlation length critical exponent \s v = 1.15(10) 
and the anomalous dimension of the magnetization is 
2^m' = 0.97(3). This latter result implies that the frac- 
tal dimension of the giant anisotropic cluster is given by: 
d'j ^ d — xjjj — 2.03(3), and evidently d'j < df. Start- 
ing with the anisotropic model we can go to the isotropic 
model by letting the couplings to be random in the ver- 
tical direction, too. Our results indicates that during 
this process the mass of the largest cluster is increasing, 
i.e. the creation of new connected parts is more effective 
than the creation of isolated sites. This result is in ac- 
cordance with the form of the phase-diagram in Fig. |51 
in which due to a similar process the ordered phase has 
a larger extent with increasing disorder. We note that in 
2d our numerical results showiS that dj^ = df, which let 
us to conjecture the exact values of the critical exponents 
in the random bond Potts model. Here we have argued 
that in 2d due to duality the creation and destruction 
processes play equivalent role and therefore the fractal 
dimension stays unchanged. 
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